BDF_coes;
Rho3=[1,data(7:9)];Sigma3=[data(10),0,0,0];
Rho4=[1,data(11:14)];Sigma4=[data(15),0,0,0,0];
Rho5=[1,data(16:20)];Sigma5=[data(21),0,0,0,0,0];
Rho6=[1,data(22:27)];Sigma6=[data(28),0,0,0,0,0,0];
theta = linspace(0,2*pi,1000);
x = exp(1i*theta);
hold on;
plot(polyval(Rho3,x)./polyval(Sigma3,x),'y','LineWidth', 2);
plot(polyval(Rho4,x)./polyval(Sigma4,x),'m','LineWidth', 2);
plot(polyval(Rho5,x)./polyval(Sigma5,x),'b','LineWidth', 2);
plot(polyval(Rho6,x)./polyval(Sigma6,x),'r','LineWidth', 2);
legend({'p = 3','p = 4','p = 5','p = 6'});
